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10  AStTRACT  (Contimf  «n  oldo  II  nocooomfy  Idottllly  Ar  AI«cA  n«MNA«f) 

An  analytic  method  for  the  reconstruction  of  the  profiles  of  refractive  indices  of  inhomog- 
eneous, stratified  regions  is  presented.  The  reflection  cisefficient,  r(K),  for  time-harmonic  electro- 
magnetic waves  is  represented  as  a rational  function  of  k,  the  wave  number.  One-dimensional 
inverse  scattering  theory  has  been  applied  to  obtain  analytic,  closed-form  expressions  for  the  pro- 
file functions,  qfxi,  from  Ka).  The  profile  reconstruction  method  is  demonstrated  with  a three- 

pole  Ka),  the  yields  q(x)  which  resembles  an  ionospheric  electron  density  profile  previously  ;vl> 
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analyzed  by  direct  scattering  methods.  The  present  communication  generalizes  previous  results  to 
oblique  incidence  and  compares  several  q(x)  obtained  from  different  rational  approximations  to  r(K). 
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INVERSE  SCATTERING  THEORY  AND 


PROFILE  RECONSTRUCTION* 

The  customary  procedure  for  constructing  a theory  for  an  electromagnetic  scattering 
phenomenon  is  first  to  assume  some  specific  model  for  the  scattering  ob>!?ct  and  then  to  calcu- 
late the  resultant  scattered  fields.  This  procedure  is  known  as  f/trinv  scattering  theory.  The  scat- 
tered fields  thus  predicted  are  compared  with  the  experimental  scattering  data  and  the  specific 
model  is  altered  until  theory  and  experiment  agree  according  to  some  acceptable  criterion. 

A different  approach  is  to  assume  only  the  general  physical  properties  of  the  scattering 
object  and  then  to  determine  analytically  the  specific  model  using  only  the  knowledge  of  the 
incident  fields  and  the  scattering  data.  This  priKedure  inverts  the  customary  analysis  of  the 
cause-and-effect  relationship  and  so  is  known  as  mwrst  scattering  theory.  This  general  problem 
can  be  simplified  if  the  scattering  object  is  assumed  to  be  an  inhomogeneous  region,  whose 
index  of  refraction  has  only  a one-dimensional  spatial  variation.  This  procedure  is  known  as 
profile  reconsiruaion. 

As  with  direct  theory,  various  approximate  methods  have  been  us.d  to  solve  inverse 
scattering  problems,  for  example  the  physical-optics  approximation*'  and  linearization 
methods*’.  Several  inverse  scattering  theories  have  been  discussed  in  recent  surveys.’ 

Miinuvripi  vihmiHcil  rcbruiir\  l.V 
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Numerical  solution  to  the  time-domain  integral  equation  for  this  problem  has  also  been  con- 
sidered/ In  the  present  communication  we  discuss  an  inverse  scattering  theory  that  may  be 
called  exact  in  the  sense  that  it  reconstructs  refractive  index  profiles  exactly  starting  from  an 
analytic  expression  for  the  reflection  coefficient. 


The  general  physical  model  which  we  consider  is  the  scattering  of  electromagnetic  waves 
from  a stratified  ionized  region,  as  shown  in  Figure  1;  this  model  has  been  used  to  study  ionos- 
pheric radio  wave  propagation.’  The  effects  of  electron  collisions  and  static  magnetic  fields 
have  been  neglected  so  that  the  relative  permittivity  of  the  inhomogeneous  region  is 


«(A:.x)  . 1 , , ^ „ 

1 r q(x),  0 


(1) 


where  t„  permittivity  of  free  space  and  k - <u/c,  the  wave  number  in  free  space,  and  where  w 
- radian  frequency  and  c - velocity  of  light.  The  profile  function,  </U),  is  proportional  to  the 
electron  density. 


The  time-harmonic  amplitude,  u{k.x),  of  the  horizontally  polarized  electromagnetic  field 
£-  satisfies  the  differential  equation  in  the  variable  .x 

u(k,x)  + - (/(jxiIwIk.jx)  - 0.  (2) 

dx^  ' ' 

where  the  spectral  variable  is  the  wavenumber  in  the  .v-direction,  k — kcos9,  u(k.x)  is  the 
Fourier  transform  of  £,.{cf.jt),  and  t — time  variables.  The  case  of  a vertically  polarized  field 
will  be  discussed  in  Section  3.4. 


The  profile  function  q{x)  in  Eq.  (2)  will  be  assumed  to  be  real,  bounded,  and  piecewise 
continuous  in  0 < .x  <oo  with  ^(.x)  s 0 for  .x  < 0.  Thus  it  is  possible  to  obtain  a solution  of 
Eq.  (2)  which  satisfies  the  asymptotic  conditions 


u{k,x) 


/•(k)p'"* e'"',  x—'—oo, 
T{K)e"‘',  .x—- (-00. 


(4) 
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The  reflection  coefficient.  rM,  and  transmission  coefficient,  7T*t),  are  assumed  to  be 
represented  by  analytic  functions  of  k in  the  complex  K-plane.  From  the  conservation  of 
energy, 

|/■(K)|^+|7'(K)|2-  1.  (5) 

where 

|r(K)P“  rM  rix)  - r(>t)  r(-  k). 

|r(K)|^-  r(K)  fuT-  hk)  n-K) 
for  real  k . and  where  /TkT  means  complex  corrugate  of  rM. 

The  inverse  scattering  problem  can  now  be  stated:  Given  the  reflection  coefficient  Hk  ) 
as  an  analytic  function  of  the  wave  number  k,  find  the  self-consistent  profile  function  ^(x). 
(Here  self-consistent  means  that  no  further  information  is  needed  to  find  a unique  profile.)  We 
will  direct  our  attention  to  the  analytic  relationship  between  the  reflection  coefficient  and  the 
profile  function.  The  problem  of  the  appropriate  data  processing  to  obtain  the  reflection 
coefficient  /-(k)  warrants  a separate  investigation 


The  profile  reconstruction  method  which  we  discuss  here  is  an  application  of  the  inverse 
scattering  theory  for  reflection  coefficients  which  are  rational  functions  of  the  wave  number." 
This  theory  was  developed  from  the  solution  to  the  inverse  Sturm-Liouville  problem.’  The 
theory  has  also  found  applications  in  quantum  scattering’  * and  non-linear  wave  propagation." 


We  will  demonstrate  the  profile-reconstruction  procedure  by  using  a third-order  rational 
approximation'  to  Hx).  The  frequency  variation  of  the  reflection  coefficient  and  the  form  of 
the  reconstructed  profile  resemble  those  that  have  previously  been  analyzed  in  connection  with 
the  solution  to  the  direct  ionospheric  scattering  problem.’  In  addition  two  more  examples  will 
be  analyzed;  r{x)  has  three  poles,  one  of  which  represents  a ”bound  state”  which  leads  to  a 
negative  profile,  and  r{x)  has  two  poles  and  one  zero,  which  leads  to  an  oscillatory  profile. 


. ■ ..1^  4., 
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These  three  examples  and  two  previously  published  examples  will  be  compared  in  order  to 
illustrate  the  relationships  between  the  general  forms  of  g{x)  and  the  pole  configurations  of 
rU), 

2 Inverse  Scattering  Problem 

The  time-harmonic  wave  amplitude  in  the  free-space  region  is 

x)  - e“"  + rMe'"".  x > 0 (6) 

where  only  the  dependence  on  the  x spatial  coordinate  is  shown.  The  corresponding  time- 
dependent  electric  field  is 

£,„(.x.cJ)  — 8(x  — ci)  -t-  R(x  + ci),  X < 0.  (7) 

where  6(x-ct)  - incident  6-function  impulse,  and  /J(x-fcf)  - reflected  transient. 

The  time-dependent  field  in  the  inhomogeneous  region  satisfies  the  differential  equation 

1 d-£,  ^ ^ 

^ - -7  ^ - qME,  - 0,  X > 0.  (8) 

ox  c'  Of 

The  retarded  electric  field  can  be  represented  in  terms  of  the  electric  field  defined  by  Eq.  (7) 
with  the  transformation"*  " 

£;  (x.cf)  - lx.cf)  ■(-  f'  K(x.()E^,  (9) 

where  the  function  K{x.  f)  also  satisfies  the  differential  equation  (8)  with  the  boundary  condi- 
tions 


0 

1 

1 

(10) 

JKix.x)  1 , 4 

* - 2 

(H) 

4 
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For  a wave  movin»  toward  ihe  right,  the  retarded  field  salisfies  the  condition 

fc',  (v.cr)  - 0.  a < 

so  that  equation  (91  together  with  equation  (7)  provides  the  integral  equation 

R ( v + cf)  -f  A ( v.cf)  + J ' A ( V.  0 « (f  + cr)  Ji  - 0 (|3) 

If  this  equation  can  be  solved  for  A(\.  f),  then  Fq  (ID  gives  the  solution  to  the  inverse 
scattering  problem 

The  integral  equation  (13)  can  be  solved  exactly  when  r(K)  is  a rational  function  of 
considered  as  a complex  variable  In  this  case  the  reflected  field  is 

A(\+i7)--^  (14) 

~ .-I 

where  Ihe  integral  represents  the  continuous  spectrum  of  r(«)  and  the  discrete  spectrum  is 
represented  by  the  sum  over  the  poles,  «„  if  any,  on  the  positive  imaginary  axis  with  the  resi- 
dues, r.  If  v(v)  > 0,  as  in  the  present  model  of  ionospheric  scattering,  then  r(K)  has  no 
poles  on  the  positive  imaginary  axis.  If  the  poles  of  /•(«)  lie  on  the  unit  circle,  then  r(<t)  is  Ihe 
mh-order  Buiierworth  approximation  and  the  integral  equation  (13)  can  be  solved'"  rather 
easily. 

The  profile  reconstruction  method  can  be  interpreted  in  terms  of  the  space-time  represen- 
tation shown  in  Figure  2.  A delta-function  impulse  is  incident  on  the  inhomogeneous  region 
and  a reflected  field  and  a transmitted  field  are  prinluced  I'sing  causality  conditions,  the  fields 
in  the  left  and  right  half-planes  can  be  reconstructed.  In  the  left  half-plane  the  total  field  is  the 
sum  of  the  incident  and  reflected  fields.  In  .he  right  half-plane  the  field  can  be  represented  by 
a linear  transformation  of  the  free-space  field  Causality  for  the  backward  light  cone  (time-like 
region)  says  that  the  reflected  field  is  not  produced  until  the  incident  field  arrives  at  v - cr,  i e 
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R (,v  + ct)  - 0 for  .V  + cl  < 0.  In  the  forward  light  cone  (space-like  region)  the  field  is  modified 
only  by  the  medium  it  has  traversed,  i.e  , A'(.v.t7)-0  for  ct  > x.  Substitution  of  this 
representation  of  the  field  in  the  original  wave  equation  provides  the  boundary  condition  which 
gives  the  profile  function 


3 Profile  Reconstiuction 


The  reconstruction  method  can  be  demonstrated  with  a third-order  rational  approximation 
to  r(K).  I e.  r(»()  has  three  poles  in  the  complex  K-plane.  The  resultant  profile  function  qix) 
resembles  electron  density  profiles  that  have  been  analyzed  by  direct  methods ^ An  alternate 
method,  which  may  be  more  convenient  if  tIk)  has  only  a few  poles,  will  also  be  demonstrated. 


3.1  Three-pole  Reflection  Coefficient 


We  consider  the  reflection  coefficient 


rix)  - 


K|*CjK, 


(k  - Kj)  (*C  — K2)  (k  — K,) 


where  kj  - -i<|  - c,  - ic;  and  k,  - -la.  The  normalization  has  been  chosen  so  that  r(0)  — 
-1.  Conservation  of  energy  requires  that 

l'•(»<)|•^l;  (16) 

this  defines  the  regions  for  the  allowed  pole  locations,  shown  in  Figure  3 by  the  shaded  portion. 
The  reflected  energy  density  |r(K)P  is  shown  in  Figure  4 as  a function  of  k for  the 
configuration  of  poles  shown  in  Figure  3;  o - 1.0,  c,  - 0.50,  Cj  - 0.499.  This  configuration 
is  also  shown  in  Figure  6.3. 


The  integral  equation  (13)  can  be  reformulated'*  without  the  assumption  ct  < x.  The 


"entire"  electric  field. 


&{x.  ct)-K(x,  ct)  + 6{x.  ct)  -E^  ix.  ct). 


NKl  Ml  MOK  ANIU'M  Kt  IHIKI  mi 


SAiities  the  integral  eviuaiiun 

Af(  v + 4/)  + »•(».  4/)  + f 4^.1 4,  O R Ct)ii(  - 0 

•'mgxi  i »*> 

and  has  ihe  properties 

f' ( v.4-f)  - - 4St\.  47).  ct  '>  \. 


(18) 


(Id) 


it\ 


(20) 


so  that  e4)uation  (20)  provides  the  solution  to  the  profile  reconstruction  problem 


By  using  the  complex  variable  s ~ <«,  integral  eituation  (15)  is  amenable  to  solution  by  a 
Laplace  transform  technique'  '•  which  will  be  used  in  this  section  The  Laplace  transform  of 
Eg  (15)  is 


where 


" f /((s)  F(\.  - 

-s)  4 A'(v,  s)  4 (i'(v.  s)  - 0. 

(21) 

•(  (s)  = r(/s)  - 

r~  R(\  + {)  c'<  Jt. 

•'O 

(22) 

A ( V.  s)  - f 

' A.  (A.  {)  c '< 

■ % 

(2.1) 

G(\.  s)  “ r 

" A ( v.  0 e '<  iii 

(24) 

In  the  present  example 


_ Aiil 

P(s)  (s  + (j)  (s'  + 24‘j  + C'"*) 
where  C'  - c,'  + c/.  In  general  we  can  write 


.V(s)  “ X " 


where  the  coefficients  » and  J,  are  real 


(25) 


(2b) 


Equation  (21)  can  be  written 

iV(s)  e"  + S(i)  F(-s.  x)  + />(s)  FU.  v)  + P(s)  0(s.  v)  - 0,  (27) 


7 


JORDAN  AND  AHN 


However  G(i.  .r)  D(s)  — /r(5.  x)  is  an  entire  function'*’  and  we  can  write,  using  (24)  and 
since  D(s)  is  a third-order  polynominal. 


where 


IX  ^ I \ ^ o ^ \ K ‘ 

e G(s.  v)  — — — - + — as  s — 

s s‘  s’ 


■ - X I-  dK(x.  ()  1 

K^-K(x.  x).  A, ^7"^ 

..  d'K(x.  0 , 

rff-’  It-- 


In  general  we  find  that 


M, 


jf  (s 


. -v)  - ^ s 1 - V 


(•>1  .-I 


where 


We  finally  obtain 


where 


A.  - 


(fKU  0 
dV 


F(x.  s)  — 


FyC.V.  s) 
Fptx.  s) 


(28) 


..2J) 


(30) 


fy(.v.  s)  — e"  ;V(s)  g(—s,  x)  — iV(s)  D(—s)  <’"''.V(i)  A^(— s)  —D(-s)  g(s,  .v). 

Fpix.  s)  - D(s)  D(-s)  - A(s)  .V(-s). 

Since  F(x.  s)  is  an  entire  function,. 


A„  - A (.v.  .v)  - lim  |s|  f(.v.  ,t)  c*'''. 
so  that  q(x)  can  be  found  from  Eq.  (11). 


(31) 


In  the  present  example 


g(s.  .v)  - -s-A„  - sy^  - y„ 


(32) 


where 
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- A^2  + (a  + 2cj)  K^  + {C^+  2cja)  K^,  (33) 

■yi  — + (fl  + 2cj)A^„. 

Using  (32)  and  (30)  in  (23),  K(x,  f)  is  found  by  an  inverse  Laplace  (ransform.  We  can  then 
solve  for  /f , and  to  obtain,  after  lengthy  but  straightforward  calculation; 

(~T^+(dj-  K„)  T + g - d,) 

K(x.  0 - (34) 


4t2(t^-  8j) 


where  the  r’s  are  solutions  of 


(t/-6j)^- 8,^-0,  7 - 1,  2.  3.  4.  (35) 

Tj  - T,.  T4  - -T,. 

8|^  “ T (o^  + 4f/)  (fl^  - 4c^). 

4 

8j  - i flJ  - (c,J  - f/). 

The  graph  of  q(.x)  for  the  rfx)  of  Fig.  3 is  shown  in  Figure  5 and  is  compared  with  other 
examples  in  Figure  6.3.  The  complete  formula  for  q(x)  has  been  previously  displayed;'  the 
analysis  leading  to  this  formula  is  summarized  here. 


3.2  Three-pole  Reflection  Coefficient  with  One  "Bound  State" 


An  example  of  a discrete  as  well  as  a continuous  spectrum  is  furnished  by  a reflection 
coefficient  with  the  pole  configuration  shown  in  Figure  6.S.  The  symmetric  poles  on  the  unit 
circle  in  the  lower  half-plane  correspond  to  the  two  symmetric  poles  for  the  third-order  Butter- 
worth  approximation,  these  poles  represent  the  continuous  part  of  the  spectral  function  for  the 
differential  equation  (2);  the  pole  on  the  positive  imaginary  axis  represents  the  discrete  part  of 


the  spectral  function: 


rM  - 


K,  - y (-s/J  - /)  , 
Kj  - -K,  . 

- i. 
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The  characteristic  Tunction,  Rix),  is  found  from  (14)  to  be 


R( 


V 1 ■*•/>/?  . 1 

x)  : f ‘ t ‘ + — e', 


(36) 


6 ' 6 ■ 3 

where  the  first  two  terms  represent  the  continuous  pan  of  the  spectrum  and  the  last  term 


represents  the  discrete  part  of  the  spectrum. 


We  will  use  this  example  to  demonstrate  an  alternate,  but  equivalent,  technique'"  " for 
solving  the  integral  equation  (13).  It  is  possible  to  construct  a differential  operator  /(p), 

p — such  that  f ip)  R(x)  - 0.  For  a three-pole  reflection  coefficient, 
ax 

f {.p)  — p'  + /(K)  + *(j  -f  *f|)  pj  " (Kl*t)  + K|»(j  -f  XjXy)  p - (37) 

so  that  in  the  present  case,  /(p)  - p'  - 1.  The  differential  operator  is  applied  to  equation  (13) 

to  obtain 


/(p)  K(x.  >')  + A (.x.  -y)  - 0. 


(38) 


and  by  symmetry 

f{-p)  K(x,  -,v)  + R(x,  y)  - 0. 

The  boundary  conditions  on  K ( v.  y)  are 

A'(.v.  y)  Iv-  » - 0. 

A'(.v.  y)  |,_  , - R'ix)  |,_o-0. 
k (v.  y)  - R (.v)  1,^-  -1. 

Eliminating  A (.v.  -y)  between  equations  (38)  and  (39)  yields  p*  K(x.  y)  - 0,  so  that 

K(x.  y)  - (',(.v)  y’  + C'4(.v)  y^  -k-  C',(  v)  y * C*3(  v)  y'  C,(.v)  y + C'o(.v). 

From  equations  (38)  to  (42)  we  obtain 


(39) 


(40) 

(41) 

(42) 


(43) 


K (.v.  y)  - - 


8v'  + 12  4.v'  + 6 2(4v'  + b)  ' 

so  that  the  profile  function  is  found  from  equation  (17)  to  be 

, . 24v(2v’  - 3)  . ,, 

9U)  - « 


(44) 


(45) 
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As  .V  — oo  ^(.v)  “ l/v’,  whuti  is  ihc  same  asymptotic  behavior  as  the  proHle  function  which 
was  derived  from  the  second-order  Butlerworth  approximation'"  There  is  a "potential  well" 

closer  to  v - 0 with  one  "bound  slate"  or  "characlcnslic  nunlc"  with  Ihc  value  1, 

corresponding  to  the  pole  on  Ihc  positive  imaginary  axis  at  k,  - + < There  is  a simple  check 
on  the  number  A/ of  bound  states  which  was  obtained  by  Bargmann'  lor  direct  uuanlum  scatter- 


ing theory 


,\f  f I'l  U ^ Af  -I  1. 


where  1/  is  the  portion  of  the  profile  function  where  </l'^  After  integrating  by  parts 

between  the  limits  0^  x we  can  evaluate  this  integral  to  xxbtain  M ^ 3-2  In  2;  so 

that  A/  “ I.  (Positive  and  negative  values  of  the  profile  function  can  be  interpreted  physically 
in  terms  of  the  scattering  o\'  vertically  polarired  waves  by  an  inhomogeneous  dielectric  region, 
discussed  in  Section  3 4 1 


3.3  Two-pole  Reflection  I'oelflcleni  with  One  Zero 


A reflection  coefticieni  with  a zero  a a-0  is  shown  in  figure  6 4 r(Kl  also  has  the 
second-order  Butlerworth  poles'" 


»C-t-/v2a  - I 

The  reconstruction  method  yields  Ihc  profile  function. 

</  \ ( X ' 

</(x>  - ^ ‘'• 

<to(X  ' 

where 


■ - 2 - -!^l(4-^/ylsin^/^x  -t  ^/^cos  v^x) 
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(>/5+>)f'-(v^-l)«--^  ^SIPn/Jv 


Since  r(K)l,_„  —0,  there  is  a potential  well  for  small  .v.  however  there  are  no  bound  states 


3.4  Vertical  Polarization 

If  the  electromagnetic  field  is  vertically  polarized  so  that  /?  - J,  //,,  then  the  timehar- 
monic  amplitude  v(A.  x)  satisfies  the  dififerential  equation 


t/  I 1 «/v(A..v) 
dx 


dx 


+ U^- 


j A^sin^O 


v(A.  x)  - 0 


(47) 


This  can  be  expressed  in  a form  similar  to  Eq.  (2)  by  using  the  local  wave  impedance  lf'(.») 
with  the  following  transformation  of  variables 


H'(.v)  - 


«(.v)  • 

4>(A..v)  — v(A..v)  yJ W' (.v)  , 

” r,  Vmo«U>  d(  . 


g(x)  - 


( ^♦”(.v) 


2 W' (v) 


W(x)  . 
2lf(.v)  ■ 


Some  values  of  H^'(.v)/(2  W'/fv))  can  cause  a negative  ^(.v).  W (.v)  can  be  found  from  Eq, 
(30)  and 


(♦'(.v)  - 


If'(O) 


[1+mo))^  ■ 

where  F{x.  0)  - F{x.  s)  \ ,_<,and  »f'(0)  - H'{x)  | ,_o. 


(48) 


4.  Discussion 


The  general  form  of  the  profile  function  q{x)  is  related  to  the  pole-zero  configuration  of 
the  reflection  coefficient  rM.  The  continuous  spectrum  is  represented  by  poles  in  the  lower 
half-plane.  The  ’smoothness’  of  q(x)  is  determined  by  the  number  of  poles  and  zeroes  of 
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r(i»):  If  r(ic)  hu  M poles  and  no  zeroes,  then  the  M-2  derivative  of  ^(.v),  and  all  lower 
derivatives  will  be  continuous  at  x - 0.  This  means  that  if  rU)  has  one  pole,  the  correspond- 
ing qix)  will  be  discontinuous  at  x - 0.  If  r(*f)  has  two  poles,  then  ^(x)  will  be  finite  at  x - 
0 but  will  have  an  infinite  slope.  If  /^(x)  has  three  poles,  both  ^(x)  and  q'ix)  are  continuous 
at  X - 0 but  there  is  an  "angle  discontinuity".  If  r(x)  has  a zero  at  x - 0.  then  q{x)  will  have 
a potential  well  since  waves  with  small  energy  penetrate  the  medium  and  are  not  reOected 
immediately.  If  a discrete  spectrum  is  present,  it  can  be  represented  by  a pole  on  the  positive 
imaginary  axis. 

The  examples  shown  in  Figure  6 illustrate  these  properties; 

1.  The  one-pole  r(x)  is  a classic  example"  and  leads  to  the  8-function  q(x)  shown  in  Fig 
6.1. 

2.  The  two-pole  case,  shown  in  Fig.  6 2.  was  previously  analyzed"'  for  the  second-order 
Buiterworth  approximation.  As  x—  oo,  q(x)  — l/.v^. 

3.  The  three-pole  example',  shown  in  Fig.  6.3.  has  the  asymptotic  behavior  q^x)  — e*"'.  It 
is  interesting  to  note  that  as  the  third-pole  x,  - -la  -*  -/oo,  the  resultant  ^(.v)  behaves 
asymptotically  like  a ^(x)  obtained  from  a general  two-pole  r(x). 

4.  If  r(x)  has  a continuous  spectrum  with  one  zero  at  x-0.  as  shown  in  Figure  6 4.  then  an 
oscillatory  profile  function  is  obtained.  However  there  is  no  discrete  spectrum  or  "bound 
state."  even  though  q(x)  does  become  negative. 

5.  If  a discrete  spectrum  is  allowed.  r{x)  will  have  a pole  on  the  positive  imaginary  axis,  as ' 
shown  in  Fig  6..S;  the  symmetric  poles  arc  taken  here  to  lie  on  the  unit  circle.  The  resul- 
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lam  profile  function  has  a "poienlial  well".  There  will  be  a "bound  stale"  or  "characterisiic 
mode  " corresponding  to  the  pole  as  the  positive  imaginary  axis  The  asymptotic  behavior 


as  .V  — oo  resembles  that  corresponding  to  the  two-pole  Bulierworih  case. 
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LUt  of  SymkoU 


X,  y,  2 
I 

c 

01 

/ 

E,  H 
E,,  H, 
E„ 

»o 

» 

k 


« 

8 

j 

( 

R 

K 


Cartesian  space  coordinates 
time  variable 

velocity  of  light  if  free  space 

angular  frequency  for  expC/ur)  lime  variation 

electromagnetic  fleld  vectors 

>>-component  of  £.  H 

Ey  in  region  .v  0 

unit  vector  in  .^direction 

incidence  angle 

scattering  angle 

wave  number  in  free  space 

wave  number  along  jr-direction  or  spectral  variable 

reflection  coefficient 

transmission  coefficient 

permittivity  of  region  x > 0 

permittivity  of  free  space 

profile  function 

Dirac  delta  function 

/K,  complex  variable  of  integration 

variable  of  integration 

reflected  transient  field  or  spectral  function 

transformation  kernel  function 
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U 

V 

W 

M 

P 

f 

& 

Xn 

r, 

‘•i 

C2 

a 

C 

A.  F.  G.  iV,  K„ 

Si  ^11 


“ Fourier  transform  of  E, 

— Fourier  transform  of  Hy. 

— scalar  amplitude  related  to  v 

— local  wave  impedance 

— number  of  bound  states 

— differential  operator 

ax 

— function  of  p 

— "entire"  electric  field 

— poles  of  rM  in  complex  K-plane 

— residues  of  rM  are  k„ 

— ReU,) 

— Im(Ki) 

— Im(K}) 

— cj*  + cf 

— auxiliary  functions  defined  by  Eqs.  (21)-(25) 
~ auxiliary  functions  defined  by  Eqs.  (26) -(3S) 
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Fi»  I - Physical  model  for  limc  harmonic  \crsioo  of  ciccirom.isnciic  «.i\c  scjucnim  h» 
inhomogeneous  lonired  region  m » > 0 Incidence  .ingle  - i),.  .md  >c,iiiering  .mgle  - 
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REFLECTED  WAVE  IN  SPACE  TIME  REPRESENTATION 
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